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ABSTRACT 

We calculate the afterglow emission for Gamma-Ray Bursts (GRBs) going off in an extremely low density 
medium, referred to as naked bursts. Our results also apply to the case where the external medium density falls 
off sharply at some distance from the burst. The observed afterglow flux in this case originates at high latitudes, 
i.e. where the angle between the fluid velocity and the observer line of sight is greater than The observed 
peak frequency of the spectrum for naked bursts decreases with observer time as t , and the flux at the peak 
of the spectrum falls off as i~ 2 . The 2-10 keV X-ray flux from a naked burst of average fluence should be 
observable by the SWIFT satellite for time duration of about 10 3 longer than the burst variability timescale. The 
high latitude emission contributes to the early X-ray afterglow flux for any GRB, not just naked bursts, and 
can be separated from the shocked inter-stellar medium (ISM) emission by their different spectral and temporal 
properties. Measurements of the high latitude emission could be used to map the angular structure of GRB 
producing shells. 

Subject headings: gamma rays: bursts - gamma-rays: theory 



1. INTRODUCTION 

A majority of long duration GRBs (lasting 10s or more) de- 
tected by the Dutch-Italian satellite BeppoSAX have detectable 
X-ray afterglows. The afterglow properties of shorter duration 
bursts is unknown, and it is possible that these bursts go off far 
away from galactic centers, where the ISM density is low. Will 
such bursts produce observable afterglows ? 

The purpose of this paper is to show that all bursts, irre- 
spective of the ISM density, should have a detectable after- 
glow emission. During the GRB the radiation is received from 
a region of the fireball of angular size along the line of 
sight to the center of explosion. Emission from higher latitudes, 
9 > r _1 , is received over a time interval that is long compared 
to the duration of the burst and, although this radiation is rela- 
tivistically beamed away from the observer, it nevertheless has 
significant magnitude. We calculate this emission and apply it 
to GRBs going off in a very low density ISM (§2). We also 
consider a case where the density of the circum-burst medium 
drops off abruptly at some radius (§3). 

2. HIGH LATITUDE EMISSION FROM A RELATIVISTIC SHELL 

Consider a spherical shell moving with Lorentz factor T. The 
shell is shock-heated at some initial time and starts to radiate. 
The emissivity in fluid rest frame, e' v ,, is a function of shell 
radius r and we assume that it is independent of the angle 9 
relative to the observer's line of sight toward the center of the 
shell. We also assume that due to the radiative electron cooling 
and the adiabatic expansion of the shell (or some other process) 
the emissivity in the observed energy band drops to zero when 
the shell radius is r c . One can show in this case that the ob- 
served peak flux at observer time t 3> r c /T 2 decreases as t~ 2 
and that the observed peak frequency decreases as t~ l . 

A simple physical explanation for these results is the follow- 
ing. The flux per unit frequency from a relativistic source mov- 
ing at an angle 9 T^ 1 is smaller by a factor of (9T) 6 com- 
pared to the case where the source is moving directly toward 
the observer (assume flat spectrum for simplicity). Integrating 



over sources located on equal arrival time surface we find the 
flux ratio to be (6T) 4 . The emission from angle 9 arrives at a 
time which is larger than the photon arrival time from a source 
at 9 = by a factor of (9T) 2 . Thus, the observer sees the 
flux falling off as t~ 2 . The observed frequency ratio in the two 
cases, for a fixed source frequency, is (9T)~ 2 oc t~ l . A deriva- 
tion for a more general case is given below. 

The flux received at frequency v from a shell moving with 
Lorentz factor T and velocity v is given by 
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where [i — cos 9, 9 is the angle between fluid velocity and the 
line of sight to the observer, and v' = vT(\ — vfi) & e' v , are the 
frequency and emissivity in the shell rest frame at radius r and 
laboratory frame time ti ab = t + rfi. For a shell of thickness Ar 
(in the lab frame) much smaller than r, the angular integration 
in the above equation can be carried out to yield 
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where p,(r) = {t lab - t)/r(t lab ). 

Let us assume that the observed spectrum during the GRB 
phase, i.e. before the shell cools, is a power-law function of 
index (3 between the observed energy band and the peak of the 
spectrum: e' v , — c'v'P, and that, as 9 is varied, the co-moving 
frame observing frequency does not pass through any breaks. 
In this case we can rewrite equation (0) as 
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As mentioned above, we consider the case where the in- 
jection of accelerated electrons stops at a certain radius r c . 
This can happen either because the internal shock has finished 
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traversing the shell or because the density of the ISM drops by a 
large factor. In both cases the electrons undergo adiabatic cool- 
ing beyond r c which leads to a sharp fall off of e' (see §3); for 
radiative electrons e' drops off even faster. 

The integrand in equation (Q) is a rapidly increasing function 
of r, hence most of the contribution comes from r ~ r c , and 
the peak flux is given by 
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where t c — t(r c ) and T c = T(r c ). The peak frequency v p of 
the observed flux decreases with time as 



u p (t) = v p {t c ) 



t + T-t, 
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For a shell that is energized by the collision with another 
shell (i.e. internal shocks) and expands in vacuum, r is also a 
measure of the duration St of the pulse emitted by the shell. For 
t 3> t c , equations (Q) and (^) become 



f» p (t) = fv P ,c (tJ > M*) = v v-,c ( 
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where f„ c = f„ (t c ) and u PtC = v p (t c ). For power-law spec- 
tra, the observed flux at a frequency v can be calculated from 
equations & @: 
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For a GRB consisting of N pulses the high latitude afterglow 
flux is the sum of flux from each pulse given in equation (0). 
For an average peak amplitude of /„ , the afterglow flux is ap- 
proximately equal to f Up N l3 ~ 1 (t G /t) 2 ~ f} ; where t G ~ NSt is 
the duration of the GRB. 

The low energy power-law index for GRBs is in the range 
of — 1 and 2 with the peak of the distribution at ~ 0, and the 
high energy index is between —3 and —0.3 with peak at —1.2 
(Preece et al. 2000). Thus the low energy light-curve for naked 
bursts is expected to fall-off as t~ 2 , whereas the light-curve at 
high energy should decline as t~ 3 2 . For synchrotron emission 
the spectral index below the peak is j3 = 1/3 and the afterglow 
from high latitude emission is expected to decline like i~ 5 / 3 
until the synchrotron peak passes through the observing band. 

As an example, consider a naked GRB lasting for 10 sec- 
onds and with a mean flux of 10 _6 ergcm _2 s _1 and the spec- 
tral peak at a few hundred keV. At 10 3 seconds after the burst 
the spectral peak is in the 2-10 keV band and the flux is 
~ 10 _12 crgcm _2 s _1 (see Figure 1). Assuming that the op- 
tical emission is not self-absorbed, the optical flux at t = 300 s 
corresponds to R ~ 25. 

The high latitude emission should be detectable by the 
SWIFT satellite as an X-ray afterglow following short dura- 
tion GRBs, which are perhaps produced as a result of neutron 
star merger in a low density medium. Some of the early X-ray 
afterglows observed by BeppoSAX have a power-law decay of 
index of > 1.6, which could have had a contribution from the 
high latitude emission. In the few cases where an X-ray after- 
glow was continuously monitored for 10 3 — 10 4 seconds after 
the main burst (i.e. GRB 910402 - Tkachenko et al. 2000, GRB 



920723 - Burenin et al. 1999, and GRB 980923 - Giblin et 
al. 1999), the X-ray light-curve exhibited a decay significantly 
slower than what is expected from the high-latitude emission, 
implying that the emission from the external shock must have 
been dominant from very early times. 

The ratio of the observed flux from the shocked ISM gas 
and the high latitude emission (Q) depends on burst parame- 
ters and, generally, on the density of the ISM. An important 
time time scale for this comparison is the shell deacceleration 
time t da 100E^ 3 (1 - i7) 1/3 (T7nor§) _1 /3 sec; where E 52 is 
isotropic equivalent of energy in observed gamma-ray emission 
in units of 10 52 erg, T 2 = To/100, T is the initial Lorentz fac- 
tor of ejecta, and -q is the efficiency factor for converting energy 
in explosion to the observed gamma-ray emission. We consider 
tda less than or greater than the GRB duration, to, separately 
below. 

Let us first consider tda ^ tG- F° r a ISM density (no) larger 
than 10 -2 cm~ 3 and for eg ~ 10 -2 , the soft X-ray domain 
is above the cooling frequency after the GRB, and in this case 
Q is independent of no. For an electron energy index p = 2, 
Q ~ 0.3(t/ G /^) /3+1 e e (t/^G) 1 " /3 (l ~ rfir)- 1 , where v G is the 
observed peak of the GRB spectrum, v is the frequency for the 
afterglow observation, and e e is the energy fraction in electrons 
in external shock. For v G jv = 20, e e = 0.1, (3 = 0, and 
r\ = 0.1, the emission from the external shock is larger by about 
a factor of 5 at the end of the GRB. This result has a weak de- 
pendence on p. 

For tda t G , as is expected for low density ISM, 

and the observing frequency smaller than the cooling fre- 
quency, and for p = 2 & (3 = 1, we find Q w 

8eTeAiX /2 CH»GM 2 vU 2 (l ~ v) 5/i V~ 5/4 ; -io is fre- 
quency in units of 10 keV. As an example, for e e = 0.1, 
e B = 10- 2 , 7] = 0.1, v G jv = 20, v w = 1, E 52 = 1 and t da = 
100s, the two emissions are equal at the deacceleration time for 
n~4x 10~ 3 cm -3 . For t da 3> t G , v greater than the cooling 
frequency, and \3 = 1, Q = 2e e (t/i da ) 3 (^ G /^) 2 (l - fy)^ 1 . 

One should be able to separate out the contributions of the 
high latitude and shocked gas emissions by using the differ- 
ence in their spectra and light-curve slopes: the X-ray spec- 
tra for shocked low density ISM is f v oc £-3(p-i)/4j,-( p -i)/2 
whereas the high-6> spectrum is the low energy part of the GRB 
spectrum, i.e. /„ oc t^^^v 15 with — 1 < < 2. 

3. FIREBALL EXPANSION INTO ISM WITH A DENSITY 
DISCONTINUITY 

In this section we consider an external shock propagating in 
a ISM which consists of two regions of different densities. The 
model we consider consists of a fireball that shock the inter- 
stellar medium, producing a standard afterglow emission, and 
then at some radius r ad the density of the medium drops pre- 
cipitously and the shell expands adiabatically so that the ther- 
mal energy of protons, electrons and magnetic field is converted 
back to the bulk kinetic energy of the shell. We follow the shell 
evolution and synchrotron radiation starting from the time of 
the free expansion of the shell. 

The thermal Lorentz factor of particles, in an adiabatically 
expanding shell, decreases as 7^ oc V -1 / 3 , and the bulk 
Lorentz factor of the shell (T) increases with time as V 1 / 3 ; 
where V = 7r0 2 r 2 Ar is the co-moving volume of the shell, 
and Ar oc r 1 / 4 is the co-moving shell thickness. 

We consider the collimation angle 8q of the ejecta to be con- 
stant in which case the thermal Lorentz factor decreases with r 
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as r 3 / 4 , and the bulk Lorentz factor of the shell increases as 
r 3 / 4 . Therefore the evolution of the thermal Lorentz factor is 



10" 3 < n 2 /n 1 < 1CT 1 . Since the observed flux for frequency 
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where T ac i and t aa - are the bulk Lorentz factor and the observer 
time respectively at the onset of the free adiabatic expansion. 
The bulk Lorentz factor T w ^ a d/lth- 

The magnetic field strength, assuming that it is tangled, de- 



creases as 



y-2/3 



1th- 



Thus the peak synchrotron fre- 



quency, in the observer frame, scales as jf h , and the peak flux 
/„ oc jth- The flux at a frequency greater than the synchrotron 
peak v m but smaller than the cooling frequency v c is given by 
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Therefore the afterglow light-curve steepens continuously; in 
the beginning of the adiabatic expansion a = 3{p— l)/4 while 
at t = 3t a d, a = 3(3p — l)/8, assuming that v m < v < v c . 
As the light-curve slope increases with time, the flux from 
higher latitudes takes over. For a shell interacting with a uni- 



form circum-burst material at r 

tad = r ad /8T 2 ad . For t > t a . 
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therefore the high latitude emission prevents a from becoming 
larger than (p + 3)/2. 

These results apply over a limited range of t. At late times 
the non-zero density of the ISM prevents the free expansion of 
the shell and the freshly shock-heated gas contributes to the ob- 
served flux. The free expansion of the shell is terminated when 
the mass of the swept-up low density gas is ~ E/T 2 d , E being 
the energy of the adiabatic shell. Thus the radius at which the 
free adiabatic expansion is terminated is r/r aa - ~ (ni/n®) 2 ' 9 , 
where n\ and n 2 are the densities of the high and low density 
ISM respectively. The time in the observer frame when the adi- 
abatic expansion ends is 
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between v m and v c is proportional to ei P ~ 1 ' > eg~ r± "' ± , values of 
e e & es for the shocked lower density medium smaller by a 
factor of 10 could reduce the flux from the low density shock 
gas so that the high-latitude emission dominates for ~ 20 t a d, 
and during this period the power-law index of the light-curve is 
a = (p + 3)/2; for v > v c , a = (p + 4)/2. 

The optical light-curve of the afterglow of GRB 00030 1C 
fell off as ~ t^ 1 for the first three days and subsequently steep- 
ened to ~ t~ 3 (Rhoads & Fruchter 2000). From simultaneous 
optical-IR observations Rhoads & Fruchter (2000) and Sagar 
et al. (2000) have found that (3 = -0.9 ± 0.1 at t ~ 4 days. 
A possible explanation for the steep decay seen at late time 
in this afterglow is that the ISM density fell off at some ra- 
dius 1 and the subsequently observed afterglow emission arose 
at 9 ^> r~\ yielding a power-law decaying light-curve of in- 
dex a = 2 — j3 — 2.9 ± 0.1, which is consistent with the data. 

Figure 2 shows the observed i?-band data for GRB 00030 1C 
and the theoretically calculated light-curve based on the model 
described here. The transition time for light-curve steepening 
is ~ 10 t a d, which is roughly consistent with the observations. 
It has been shown by Kumar & Panaitescu (2000) that the time- 
scale for light-curve steepening due to jet edge effects in a ho- 
mogeneous ISM is roughly comparable. The late time power- 
law index according to the jet model is a = 1 — 2/3 = 2.8 ± 0.2 
for adiabatic electrons radiating at optical frequencies, which 
is also consistent with the data. The different relationship be- 
tween a and (3 in these two models can be used to distinguish 
between them. 

4. CONCLUSIONS 

The main conclusion of this work is that gamma-ray bursts 
going off in vacuum - naked GRBs - should have X-ray af- 
terglow emission detectable by the X-ray telescope aboard 
SWIFT satellite for about an hour after the GRB. This radiation 
originates at the high latitude, 9 3> part of the gamma-ray 
emission surface. The flux in a fixed observer energy band be- 
low the peak falls off as t~ 5 / 3 , while the peak flux decreases as 
t~ 2 . The peak frequency of the observed flux falls-off as t . 

For a burst going off in a non-zero density ISM the early 
afterglow flux, within the first hour, is the sum of emission 
from high latitude and the shocked ISM. The two can be dis- 
tinguished based on the differences between their spectral and 
temporal slopes; the X-ray spectrum for the shocked ISM is 
~ v^ 1 , whereas for the high latitude radiation the spectrum 
should be the same as the GRB spectrum at low energies, i.e. 
~ v - 3 . The measurement of the high latitude afterglow emis- 
sion should help map the irregularities in the ejecta producing 
the GRB and their collimation before these are detected in the 
emission from the shocked ISM. 

The radiation emitted from latitudes 9 > sets an upper 
bound on the steepness of the flux decline: we expect the ob- 
served 7-ray flux for each individual peak within the burst to 
fall off less rapidly than t~( 2 ~^\ where t is measured from the 
peak of the pulse and (3 is the spectral index (f u oc v@). A more 
rapid flux decline would be an indication of either an extremely 
small jet opening-angle or a very inhomogeneous shell, as in 
the model suggested by Kumar & Piran (2000). 

Another straightforward consequence of the high latitude 
emission is that the power-law decline for the afterglow light- 



For nz/nx =0.1 (0.01) free expansion is terminated at 
t f /tad = 2.8 (4.2). 

The value of a reverts back to 3(p — l)/4 when the emis- 
sion from shocked low-density ISM takes over. At the time 
when the adiabatic expansion of the shell ends the ratio of the 
flux from the low density shocked gas to the flux at t a d is 
~ (n 2 M) (7 ~ p)/12 - If the fractional energies in electrons, e e , 
and magnetic field, es, are same for shock in the high and low 
density ISM then the flux at tt due to high latitude emission 
and the low density shocked gas are approximately equal for 

1 The optical emission of the afterglow of GRB 000301C shows considerable variability prior to the steepening of the light-curve decline. This suggests that there 
are significant fluctuations in the ISM density, and so it is not altogether surprising that the external medium density drops to a small value at some radius. 
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curve can not be larger than about 3, even when the fireball 
expands into vacuum. The observed late time power-law index 
for the light-curve of GRB 000301C is about 3, which is larger 
by about 2 compared to the early time index. This large and 
rapid steepening of the light-curve could arise when the late 



time light-curve is dominated by emission from high latitudes. 

We thank Bohdan Paczyriski and Tsvi Piran for useful dis- 
cussions. 
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FIG. 1. — 2-10 keV light-curve arising from high latitudes, i.e. 8 > r _1 . The ISM density is assumed zero. The peak flux for the GRB is taken to be 
10 -12 erg cm _2 s _1 cV _1 ; the burst duration is 10 seconds, the low energy slope of the spectrum is 1/3 and the high energy index is —1. The burst duration and 
the peak flux set the x- and y- axis scales, and the spectral slopes set the power-law decline of the light-curve (see text). The power-law index of the light-curve is 5/3 
for t < 400s and steepens to 3 when the peak frequency passes through the x-ray band. The inset shows the GRB light curve in 20-500 keV band, the mean pulse 
width is 0.4s. 
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FIG. 2. — i?-band light-curve (continuous line) from a spherical remnant running into a uniform density ISM which ends at some radius. Beyond this radius, at 3 
day in the observer frame, the shell undergoes adiabatic expansion. Before the break the optical flux decay slope is 3(p — l)/4, and after the break it is (p + 4)/2 
for frequencies above the cooling break. The parameters for the model are: E = 2 X 10 52 erg, p = 2.6, n = 1 cm~ 3 , eg = 5 X 10~ 3 , t e = 0.1. The cooling 
frequency crosses i?-band around 5 days. The data is taken from Sagar et al. (2000), Massetti et al. (2000) and from GCN Circulars: Bernabei et al. 2000, Bhargavi 
& Cowsik 2000, Castro-Tirado et al. 2000, Frachter et al. 2000, Fynbo et al. 2000, Gal- Yam et al. 2000, Garnavich et al. 2000, Halpern et al. 2000, Mujica et al. 
2000, Veillet et al. 2000 . 



